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^ I review some recent progress and open problems in the calculation of heavy- 

■ quarkonium inclusive decay widths into light particles in the framework of QCD 

' non-relativistic effective field theories. 
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. 1. Introduction 

^ • In recent years, the study of heavy quarkonia (66, cc, ...) has gone through 

several theoretical and experimental advances. 

From the theoretical side the introduction of non-relativistic effective 
. field theories (EFTs) of QCD^'^ has put our description of these systems on 

the solid ground of QCD. It has made it possible, in the case of several ob- 
servables, to factorize the high-energy dynamics into matching coefhcients 
calculable in perturbation theory and the non-perturbative QCD dynamics 
' into few well-defined matrix elements to be fitted on the data or calculated 

on the lattice. Systematic improvements are possible, either calculating 
higher-order corrections in the coupling constant or adding higher-order 
operators. 

Heavy quarkonium, being a non-relativistic bound state, is characterized 
by a hierarchy of energy scales to, mv and mv^ , where to is the heavy-quark 
mass and v <^ \ the relative heavy-quark velocity. A hierarchy of EFTs 
may be constructed by systematically integrating out modes associated to 
these energy scales. Integrating out degrees of freedom of energy to, which 
for heavy quarks can be done perturbatively, leads to non-relativistic QCD 
(NRQCD)^. This EFT still contains the lower energy scales as dynamical 
degrees of freedom. In the last years, the problem of integrating out the re- 
maining dynamical scales of NRQCD has been addressed by several groups 
and has now reached a solid level of conceptual understanding (an extended 
list of references can be found in^^). The ultimate EFT obtained by sub- 
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sequent matchings from QCD, where only the hghtest degrees of freedom 
of energy mv"^ are left dynamical, is called potential NRQCD, pNRQCD^. 
This EFT is close to a quantum-mechanical description of the bound system 
and, therefore, as simple. It has been systematically explored in the dy- 
namical regime Aqcd < rnv^ in^'^ and in the regime mv^ <^ Aqcd ^ rnv 
jji4,6,7 rpjjg quantity Aqcd stands for the generic scale of non-perturbative 
physics. 

From the experimental side new data have recently been produced for 
heavy- quarkonium observables. Measurements relevant to the determina- 
tion of heavy-quarkonium inclusive decay widths have come from Fermilab 
(E835)^ BES^ CLEQiO'" and Bellei^. 

In the following I will review recent progress in our theoretical under- 
standing of inclusive and electromagnetic heavy-quarkonium decays. I will 
recall the NRQCD factorization formulas in Sec. 2 and the pNRQCD fac- 
torization in Sec. 3. The presented pNRQCD formulas apply to quarkonia 
that fulfil mv > Aqcd ^ mv"^. In Sec. 4, I will discuss the main difficul- 
ties that still prevent us from obtaining precise predictions from the above 
formulas and indicate what progress has been made towards an eventual 
solution of these diSiculties. 



2. NRQCD factorization 

The NRQCD factorization formulas are obtained by separating contribu- 
tions coming from degrees of freedom of energy m from those coming from 
degrees of freedom of lower energy. In the case of heavy-quarkonium decay 
widths, the first are encoded in the imaginary parts of the four-fermion 
matching coefficients, f, gi,s.ee.-/i(^^~^^ Lj) and arc ordered in powers of 
as{m). The second are encoded into the matrix elements of the four-fermion 
operators on the heavy-quarkonium states \H) {{...)h = {H\ . . .\H)). 
These are non-perturbative objects, which are counted, at least, in pow- 
ers of mv. Matrix elements of higher dimensionality are suppressed in v. 
Including up to the NRQCD four-fermion operators of dimension 8, the 
NRQCD factorization formulas for inclusive decay widths of heavy quarko- 
nia into light hadrons {LH) read^'^^: 

+lmfsCSi) {OsCSi))vcAnS)+lmfsCSo) {OsC So))vQinS) 
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r{PQ{nS) ^LH) = ^ yimfiCSo) {OiC So)) P^inS) 

+lmfsCSo) {OsCSo))p^^nS)+'i^fs{'S,) (08('5i))p«(„s) 
+lrag,CSo) +Im/sCPO ^ (2) 

nXQinJS) -^LH) = ^ (lmf,C'+'Pj) <Oi^"'^'Pj))x.inJS) 

+/8r+'5s)(08('5o)>xa(n7S)). (3) 



At the same order the electromagnetic decay widths are given by: 
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T{VQ{nS) ^ e+e-) = ( Im/eef^i) (OEM('^i))yQ(nS) 



nXQinJl) - 77) = 2ImA,(3p,) (O^M^P/p^.M!) ^ ^ ^ .(6) 

The symbols Vq and Pq indicate respectively the vector and pseudoscalar 
S-wave heavy quarkonium and the symbol xq the generic P-wave quarko- 
nium (the states xq("10) and XQ(nJl) are usually called hQ{{n — l)P) and 
XQj{{n - 1)P), respectively). 

The operators 0,Vi^s,em(^^~^^Lj) are the dimension 6 and 8 four- 
fermion operators of the NRQCD Lagrangian. They are classified in depen- 
dence of their transformation properties imdcr colour as singlets (1) and 
octets (8) and under spin (S), orbital (L) and total angular momentum (J). 
The operators with the subscript EM are the singlet operators projected on 
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the QCD vacuum. The exphcit expressions of the operators may be found 
ini'7. 

The imaginary parts of the four-fermion matching coefficients have been 

calculated over the last twenty years to different levels of precision. In the 
following I will indicate, to the best of my knowledge, some recent literature 
where their updated value may be found. If the case arises, references to the 
original literature may also be found there. The imaginary parts of fg,{^Si), 
fsi^So), fs^Pj), and fXSo) (this originally calculated in^O) 

have been calculated up to order in^^. A different result for /i('^Po) 
and /i(^P2) is in^^. The imaginary part of /i(^S'i) has been calculated up 
to order in^^, the imaginary part of gi(^S'i) at order may be found 
in^^, the imaginary part of gi{^So) at order in^ and the imaginary part 
of fsi^Pi) at order in^^. Where the electromagnetic coefficients are 
concerned, the imaginary part of fed^Si) has been calculated up to order 
a^ckg in"'^^, the imaginary parts of f^^{^So) and f^'y{^Pa.2) up to order a^as 
(originally calculated in^°'-^^) and geei^Si) and g-y-yCSo) up to order may 
be found in^. 

The NRQCD matrix elements arc poorly known. They may be fitted 
on the experimental decay data, as in^^, or calculated on the lattice^^. The 
matrix elements of singlet operators may be linked at leading order to the 
Schrodinger wave functions in the origin^ and are often evaluated by means 
of potential models^'^. 

3. pNRQCD factorization 

The pNRQCD expressions of the NRQCD matrix elements arc obtained by 
integrating out degrees of freedom of energy larger than mv^. Four different 
situations are possible: Aqcd ^ mv'^, Aqcd ~ mv'^, mv ^ Aqcd ^ mv'^ 
and Aqcd ~ mv. In the first situation the NRQCD matrix elements may 
be calculated perturbatively up to non-perturbative corrections of the form 
of local condensates. This situation may apply to the ground state of bot- 
tomonium. Explicit formulas have been worked out for the electromagnetic 
decay of the T(15) in^^. In the situation Aqcd ^ mv'^ the NRQCD matrix 
elements turn out to be convolutions of non-local condensates and pertur- 
bative Green functions. In this situation, which may apply to the ground 
states of bottomonium and charmonium, only the real part of the spec- 
trum has been studied'''. Calculations of decay widths have not been done 
so far. In this situation, however, the state and flavour-dependent part 
will not factor out from the integrals. The situations mv ^ Aqcd ^ mv'^ 
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and Aqcd ^ "mv have been studied in^"'''^. In principle they may apply to 
all heavy quarkoniuin states below threshold and above the ground state. 
Both situations give the same physical results and will be summarized in 
the following. The main point here is that, since Aqcd is well separated 
from the energy scale of the bound state, mv"^ , the state dependence of the 
NRQCD matrix elements factorizes in the quarkonium wave function, while 
contributions coming from excitations of order Aqcd are encoded into a few 
universal constants. 



3.1. pNRQCD factorization in the case mv > Aqcd^tuv^ 

The matrix element of the NRQCD operator O on the heavy-quarkonium 
state \H) at rest, P = 0, with quantum numbers n, j, I and s, may be 
written in terms of the eigenstates |k; X1X2) of the NRQCD Hamiltonian as^ 

iO)H = ^ o|p ^ / d'r J Sr' J d^R J d^R' (P = 0|R)(ry7s|r) 



(0;xiX2| jd^^O{i) 



(R'|P = Q){v'\njls) , 



where r = xi — X2, r' = x'^ — Xj, R = (xi + X2)/2 and R' = (x'^ + X2)/2. 

In the situation mv > Aqcd ^ rnv"^ and assuming that all the higher 
gluonic excitations of the two heavy quarks develop a mass gap of or- 
der Aqcd, the matrix element (0;xiX2|O(^) |0;x']^X2) has been calculated 
order by order in 1/m in^^'^. Once normalized to m, we get up to 
X {Aqijj^ / 771"^, E/m)) for the NRQCD S-wave matrix elements and 
up to 0{v^) for the NRQCD P-wave matrix elements^: 

(OiiSi))v,inS)=CA^^\l-—^ + ^ + ^j, (7) 

1^(1^ .. _r Ko(0)l' (, E^o 2^3 , 2f i^'*) 

{Ui( bo))pQ{nS) - (^A ^ |1- 

/n f3Q\\ n l^^noMl^ M -^no 2^3 

(Cem( -^ij/VQCnS) = I 1 



in -r l-^no(Q)l^ I 1 2^3 , 2g;^^'^^^ ^pL 

{Oem{ i'o))pg(„s) - Ca — — — I 1 - — — + ^_2 + 




''Possible threshold effects have been neglected. As suggested in-^^ they may be large for 
the V'(25). 
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(10) 

{OlC'^'Pj)).,inJS) = (Oem(^^+^P.)>,,M5) = l^lR^I'm', (11) 

^ TT 

{'PlCSi))vQ{nS) = (^i(^S'o))pq(„S) = {'PEMCSi))vQ{nS) 

= (Pem('5o))p^(„5) = cj-^^^ (mEl^^ - £,) , (12) 

(08('''S'i))vq(„S) = {08CSo))pQ{nS) 

- 27r 3^;^^ ' ^^^> 



/naQW (08(^'S'i))po(nS) 

(Csi J0)}VQ{nS) = ^ 



(08('P7)>y«(nS) 



|i?i°J(0)P / (CV2-CV)c|Si' 
= ^^^^ [ 3^^ ) ' ('^^ 

{08(^Pl))pQ(nS) 

= (15) 



{Us{ ^0))xQ{nJS) - ^ —2 ^3, 



(16) 



where R^q = R^q{1 + 0{v)) is the radial part of the vector S'-wave function, 
i?^Q = -R^o''(l + 0{v)) the radial part of the pseudoscalar S-wave function, 
r!^^ ' the derivative of the LO P-wave function and E^^q ~ MH — 2m ^ mv"^ 
the LO binding energy. A consistent way to calculate the wave functions is 
from the real part of the pNRQCD Hamiltonian obtained up to order l/m^ 
in^. Any other 5- wave dimension-6 matrix element is at NNLO and any 
other S'-wave dimension-8 matrix clement is at LO. The symbol cp indi- 
cates the NRQCD chromomagnetic matching coefficient, which is known at 
NLO from^''. All the non-perturbative dynamics is encoded into the wave 
functions and the six chromoelectric and chromomagnetic correlators Ei, 
£3, £^^\ £'i''^\ £'^'^^^ and Bi. The exact definition of the correlators may 
be found in^. It is important to note that, owing to the factorization of the 
state and flavour dependence into the wave function, on the whole set of 
heavy-quarkonium states below threshold the number of non-perturbative 
parameters needed to describe inclusive decays has diminished with respect 
to NRQCD, so that definite new predictions are possible. In practice, one 
may consider ratios where the wave- function dependence drops out^^''' and 
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fix the correlators either through lattice calculations or through specific 
models of the QCD vacuum^* or on some set of experimental data. In 
particular, new predictions for P-wave bottomonium inclusive decay- width 
ratios were made in this way in^^'^^ before the CLEO-III data were made 
available^^. 

4. How to use the above formulas 

Here I would like to point out that there are two potential problems that 
may still prevent us from getting precise determinations of heavy-quarko- 
nium decay widths from the above formulas. These problems seem to be 
often underestimated in the existing literature. Eventually, they may lead 
to an underestimate of the errors associated to the calculated values of the 
decay widths or to the extracted values of the matrix elements and Q:s(m). 

i) It has been discussed, in particular in^° (but see also"^^), that higher- 
order operators, not considered in the above formulas, can be numerically 
quite relevant. This may be the case particularly for charmonium, where 

^ 0.3, so that relativistic corrections are large, and for P-wave decays 
where the above formulas provide, indeed, only the leading-order contribu- 
tion in the velocity expansion. In order to be specific, in the case of the Xc 
decay width into two photons the corrections of relative ©(Aggj-j/m^, E/m) 
to the formulas (6) are given by^°: 

rric \ I 

1.0^_4Bf^V 

rric \ rn^ I 

where ^ (Aqqj-,/to^) and ;B (Aq^q/to^) are some specific combinations 
of matrix elements of order AQ^j-j/m^. Choosing E^^/rric — 0.3 and 

B(Aqqj3/toJ?)| ~ 0.1 (or 0.3)'' the relative corrections may be — (0.7± 
0.3) (or -(0.7±0.9)) for r(xcO ^ 77) and -(0.3±0.4) (or -(0.3±1.2)) for 
r(Xc2 77) and, therefore, potentially as large as the leading one. The 
two ± values refer to the two possible choices of sign for A and B. Note 
that for some specific choice of sign, large cancellations may also occur. 

The introduction of higher-order matrix elements may spoil the predic- 
tive power of the NRQCD factorization. A way out is provided by pN- 



'^r(xco - 


* 77) 


r(o)(xco- 


-^77) 


sr{xc2 - 


■> 77) 


rw(xc2 - 


77) 



''The first choice corresponds to a hierarchy of the type Aqcd ~ mv^, the second of the 
type Aqcd ~ mv. 
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RQCD, which can be used to reduce also these new matrix elements to a 
few universal correlators'". 

ii) The convergence of the perturbative series of the four-fermion match- 
ing coefficients is often bad. Let us consider, for instance, the following 
matching coefficients (uf = 3): 



ImfsCSo) = 
ImfsCSi) = 
Im/i(3Po) = 
Im/i(3p2) = 



...)x 



...)x 



eta 

1 + 11.1 — 

IT 

1 + 13.7^), 
1 + 10.3^V 

TT / 

1 + f 13.6 - 0.44 log 



1 - 



(o.73+ 1.67 log— ^ 



2m J 

JL 
2m J 



TT / 



Apart from the case of Im/i(3p2), the series in of the other coefficients 
does not seem to converge. This behaviour may not be adjusted by a suit- 
able choice of the factorization scale /z, which enters only in Im/i(^Po); in 
this case lowering the factorization scale below 2m would make the con- 
vergence worse.'' A solution may be provided by the resummation of the 
large contributions in the perturbative series coming from bubble-chain di- 
agrams. This analysis has been successfully carried out in some specific 
cases in'^^. A general treatment is still missing, in particular in the case of 
P-wave decays. 

Finally, I would like to note that neither contributions coming from 
higher-order matrix elements of the type described in (i) nor resummations 
of large contributions in the perturbative scries of the matching coefficients 
discussed in [ii) have been considered in recent determinations of (e.g. 
iji32,2i^ and of the NRQCD matrix elements (e.g. in^^) from the charmo- 
nium P-wave decay data. 



5. Conclusions 

The progress in our understanding of non-relativistic effective field theories 
makes it possible to move beyond ad hoc phenomenological models and have 
a unified description of the different heavy-quarkonium observables, so that 
the same quantities determined from a set of data may be used in order to 

'^In the case of the coefficient Im/ee Si ) , which presents a similaj: problem, it has been 
noted by Beneke et al. in^^ that in the bottomonium case a suitable choice of the 
fax;torization scale may adjust the convergence of the series up to order a^a^ . 
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describe other sets. Moreover, predictions based on non-relativistic EFTs 
are conceptually solid, and systematically improvable. Therefore, the re- 
cent progress in the measurement of several heavy-quarkonium observables 
makes it meaningful to ask whether their precise theoretical determination 
is feasible. 

Here, I have focused on inclusive and electromagnetic heavy-quarkonium 
decays. In the framework of NRQCD, heavy-quarkonium decay widths may 
be expressed in terms of matrix elements, which, in the case of matrix ele- 
ments of singlet operators, may be linked to the wave function in the origin. 
This is also the case of some matrix elements entering in the description 
of production processes^. In the framework of pNRQCD also some octet 
matrix elements may be expressed in terms of the wave function in the 
origin and some universal non-local gluon-field correlators. These correla- 
tors enter also in the expression of the masses of some heavy-quarkonium 
states^'"^^. Specific predictions have been discussed in^^^^. There are, how- 
ever, still some difficulties that hinder in several cases precise predictions 
of heavy-quarkonium decay observables. These are essentially related to 
the control of higher-order corrections in the velocity and ctg expansion and 
have been addressed in the last section. In principle, the tools to overcome 
these difficulties already exist, so that progress in the field is expected from 
the coordinated effort of the heavy-quarkonium community in the near 
future^^. 
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